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' Abstract 

^ ■ In a recent string theory motivated paper, Nicohni, Smailagic and Spallucci (NSS) presented 

0^ ■ an interesting model for a noncommutative inspired, Schwarzschild-hke black hole solution in 

, 4-dimensions. The essential effect of having noncommutative co-ordinates in this approach is 

' to smear out matter distributions on a scale associated with the turn-on of noncommutativity 

' which was taken to be near the 4-d Planck mass. In particular, NSS took this smearing to be 

o ■ 

. essentially Gaussian. This energy scale is sufficiently large that in 4-d such effects may remain 

o 

invisible indefinitely. Extra dimensional models which attempt to address the gauge hierarchy 
problem, however, allow for the possibility that the effective fundamental scale may not be far 



* t 



Oh 
p. 



1^ ■ from ~ 1 TeV, an energy regime that will soon be probed by experiments at both the LHC and 



X 



ILC. In this paper we generalize the NSS model to the case where flat, toroidally compactified 
extra dimensions are accessible at the Terascale and examine the resulting modifications in 
black hole properties due to the existence of noncommutativity. We show that while many of 
the noncommutativity-induced black hole features found in 4-d by NSS persist, in some cases 
there can be significant modifications due the presence of extra dimensions. We also demonstrate 
that the essential features of this approach are not particularly sensitive to the Gaussian nature 
of the smearing employed by NSS. 



*Work supported in part by the Department of Energy, Contract DE-AC02-76SF00515 
^e-mail; "rizzo@slac.stanford.edu 



1 Introduction and Background 



The theoretical effort that has gone into understanding the fuh details of string/M theory has 
inspired a number of ideas which, on their own, have had a significant impact on particle physics 
model building and phenomenology. One of the more recent developments of this kind has been 
the resurgence^r of interest in noncommutative (NC) quantum field theories[2] and, in particular, 
the question of how a NC version of the Standard Model (SM) may be constructed and probed 
experimentally [1] . 

The essential idea behind NC constructions is that the commutator of two spacetime coor- 
dinates, now thought of as operators, is no longer zero. In its simplest form, for a space with an 
arbitrary number of dimensions, D, this is can be written explicitly as 

[xa-,xb\ = iOab = ^-j^ ' (1) 

where A^vc is the mass scale associated with NC and cab is normally taken to be a frame- 
independent, dimensionless anti-symmetric matrix with constant, real, typically 0(1) entries; it 
is not a tensor. Here we assume that vastly different NC scales do not exist depending upon the 
values oi A,B. ^ In a general string theory context one might imagine that A^c would naturally 
not be far from the 4-d Planck scale, Mpi, and that the cab are generated due to the presence 
of background 'electric' or 'magnetic' type fields. Most of the phenomenological studies of NC 
models ^,4 have assumed that we live in 4-d and that A^vc ~ 1-10 TeV so that we have access to 
this scale at, e.g., the LHC or ILC. However, if the NC scale A^c is indeed of order Mpi, then 
probing NC physics directly may prove difficult in the near term. 

One way of possibly observing NC is its effects on the properties of black holes (BH). 
In order to analyze this problem at a truly fundamental level one would need to successfully 

construct the NC equivalent of General Relativity. Attempts along these lines have been made 

*In what follows, upper case Roman letters will label all D dimensions while Greek (lower case Roman) letters 
will cover the range 0-3 (5 to n + 4) . 
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in the literature|S] but no complete and fully compelling theory of this type yet exists. Recently, 
Nicolini, Smailagic and Spallucci (NSS) [S] have considered a physically motivated and tractable 
model of the possible NC modifications to Schwarzschild BH solutions. The essential ideas of this 
picture are: {i) General Relativity in its usual commutative form as described by the Einstein- 
Hilbert action remains applicable. This seems justifiable, at least to a good approximation, if NC 
effects can be treated perturbatively. The authors in Ref.i5 have indeed shown that the leading 
NC corrections to the form of the Einstein-Hilbert action are at least second order in the 9ab 
parameters, (ii) NC leads to a 'smearing' of matter distributions on length scales of order ~ ^jvc" 
Thus the usual '(5-function' matter source of the conventional Schwarzschild solution is replaced by 
a centrally peaked, spherically symmetric (and time-independent) mass distribution which has a 
size of order ~ ^nxj- This, too, seems justifiable based on the results presented in Ref.fH], which 
note that matter actions from which the stress-energy tensors are derived are modified at leading 
order in the Bab parameters. Based on earlier workjjj, NSS took this smeared distribution to be 
in the form of a spherical 3-d Gaussian in 4-d whose size, due to the spherical symmetry, was 
set by a single parameter, 6, indicative of the NC scale. Though such a picture leads to many 
interesting properties for the resulting BH (to be elaborated on below), since Mpi ~ A^vc was 
assumed, as would be natural in 4-d, such BH are not immediately accessible to experiment or to 
direct observation. 

Another interesting prediction of string theory is that several extra dimensions must exist. 
However, extra spatial dimensions, in models with an effective fundamental scale now in the 
TeV range, have been discussed as possible solutions to the hierarchy problem |H1 1^1. In the case of 
'flat' extra dimensions, e.g., in the model of Arkani-Hamed, Dimopoulos and Dvali (ADD)|S], the 
4-d Planck and fundamental scales are related by the volume of the compactified extra dimensions: 

Mil = VnM:+^ , (2) 
where Vn is the volume of the compactified manifold. Assuming for simplicity that these extra 
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dimensions form an n-dimensional torus, if all compactification radii {Rc) are the same, then 
Vn = (27r i?c)"- In such a scenario gravity becomes strong at M,, and not at Mpi which is viewed 
as an artifact of our inability to probe gravity at scales smaller than Rc. This scenario has gotten 
a lot of attention over the last few years and the collider phenomenology of these types of models 
has been shown to be particularly rich 10 . In such a scheme it would be natural that the NC scale, 
^NCi would now also be of order ~ TeV allowing it to be accessible to colliders. Furthermore, 
the copious production of TeV-scale BH at colliders also becomes Dossible|ll| and the nature of 
such BH could then be examined experimentally in some detail. Thus it is reasonable to ask if 
the properties of such TeV-scale BH may be influenced by NC effects, which originate at a similar 
scale, and if these effects are large enough to be observable in collider data. 

The goal of the present paper is to begin to address these issues. In particular we will 
examine how NC BH in D dimensions differ from those in 4-d as well as from the more conventional 
L'-dimensional commutative Schwarzschild BH traditionally analyzed at colliders. Furthermore, 
we will demonstrate that the essential features of this scenario are not particularly sensitive to the 
detailed nature of the NC smearing. 

2 Analysis and Results 

We begin our analysis by reminding the reader that we will assume that D = 4 + n-dimensional 
gravity, and BH in particular, can still be described by the conventional Einstein-Hilbert (EH) 
action, i.e., 



with R being the Ricci scalar and being the (reduced) fundamental scale as appearing in the 
ADD relationship above. It is important to recall that for the ADD scenario with n > 2, M* is 
only weakly constrained by current collider experiments, i.e., one finds that M* > 0.38 — 0.60 TeV, 
depending on the value of n, when the bound on the more commonly used GRW^U] parameter 




(3) 
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Mu > 1.5 TeV is employed. 

In the present and NSS approaches the basic effect of NC is proposed to be the smearing 
out of conventional mass distributions. Thus, following NSS0, we will take, instead of the point 
mass, M, described by a 5-function distribution, a static, spherically symmetric, Gaussian-smeared 
matter source whose NC 'size' is determined by the parameter \/6 ^ ^^c'^ 



M 



(4vr0)("+3)/2 



e-^'l^'. (4) 



Here will we explicitly assume that both the horizon size of our BH and the NC parameter ^/O are far 
smaller than the compactification scale Rc so that the BH physics is not sensitive to the finite size of 
the compactified dimensions. For ADD-type models this can be easily verified from the numerical 
results we obtain below as the BH horizon size will typically be of order ~ \/M^ while Rc is generally 
many orders of magnitude larger as long as the value of D is not too large|12j as is certainly the 

case when D < 11. Note that the value of ^fO is directly correlated with the NC scale and is 
certainly proportional to it; however, within this treatment the exact nature of this relationship 
is unspecified and would require a more detailed model to explicitly determine. It is sufficient to 
remember only that ~ 1/ K^c- It is important to realize that many such parameterizations of this 
peaked smeared mass distribution are possible which should lead to qualitatively similar physics 
results. However, the various predictions arising from these may differ only at the 0(1) level or 
less, as long as the detailed structure of the peaked mass distribution is not probed. We will discuss 
this issue further below. 

The metric of our I?-dimensional space is assumed to be given by the usual D-dimensional 
Schwarzschild form 

ds^ = e^dxl - e^'dr^ - r^dfl^^^ ■ (5) 

Here we will be searching for Schwarzschild-like, spherically symmetric and static solutions with 
1/ and fi being functions only of the co-ordinate r and we will further demand that e'^'^ ^ 1 as 
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r —>■ oo; this will require that u = — /x in the solutions of Einstein's equations as in the usual 
commutative scenario since the EH action and resulting field equations remains applicable. Note 
that the surface of the sphere, can be simply described by a set of — 2 = n + 2 angles, (pi, 

where i = 1, ..,n + 2. 

Given the assumed form of the matter density, p, above and the expectation that v = — /x 
as usual, two components of the diagonal stress-energy tensor, T^^ , are already determined, i.e., 
= Tq = p. As noted by NSS, the remaining components, T? (no summation), all tt, + 2 of which 
are identical due to the spherical symmetry, can be obtained from the requirement that T"^^ have 
a vanishing divergence, T"^^;b = 0, where the semicolon denotes covariant differentiation. Given 
the nature of our metric it is easily seen that both T°^;o = and r**;j = 0, for all i, automatically. 
The remaining equation T^^;r = then yields the explicit result 



(6) 



Since by construction = Tq and, noting that g^dj-gu = 2/r, for all (unsummed) i, this yields 



for all i (without summation) . This reproduces the 4-d NSS result in the limit when n — 0. 



(7) 



With our metric the non-zero components of the Ricci tensor are given by (with the index 
i not summed) 



R] = 



e''{l + n + riy')-{n + l) 



(8) 



where now a prime denotes partial differentiation with respect to r. The Einstein equations resulting 
from the EH action augmented with our matter distribution can be conveniently written in the form 
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where T is the trace of the stress-energy tensor, T = T^. Note that given our assumptions there 
are only two distinct Einstein equations. Writing 500 = e'^ = 1 — A{r), the i?* Einstein equation 
leads to the following first order differential equation for A{r): 

from which we obtain the following solution, after substituting the above expression for p and 
demanding that A{r) — as r — 00: 



A{r) = 



1 M 1 r'/'^ 



M*'*+2 (n + 2)7r("+3)/2 r 



^ / dt e-*t(-+i)/2 . (11) 
^0 



This result is seen to reduce to that of NSS when n — > as well as to the usual D-dimensional 
Schwarzschild solution when 6^0. The remaining Rq Einstein equation just returns to us the con- 
tinuity equation for the stress-energy tensor, hence, nothing new. Note that given our assumptions 
these results allow us to calculate A{r) for any chosen form of the mass distribution p{r). 

The horizon radius, Rh, occurs at values of r where 500 = 0; ^-e., where A{r) = 1. Defining 
for convenience the dimensionless quantities m = M/M^,, x = M^Rh and y = M^^fO, along with 
the constant 

_ (n + 2)7r("+3)/2 



r(n±3) ' 

the horizon radius can be obtained by solving the equation 



(12) 



x"+i = -F„(z), (13) 
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where z = x/y and the functions Fn{z) are given by the integrals 

= / dt e-*t("+i)/2 . (14) 

These integrals be performed analytically when n is odd, e.g., Fi = 1 — e~'^(l + q), F3 = 1 — e~'^(l + 
q + and F5 = 1 — 6^*^(1 + g + + g^/6), eic, with g = (For n even these functions 

can be expressed in terms of combinations of error functions.) Given their definition it is clear that 

the Fn vanish as ~ 2""*"'^ when z ^ and they are seen to monotonically increase with increasing 
z; as z — > cxD our normalization is such that F„ — > 1. This implies that the BH mass, m, diverges as 
either x ^ 0, 00 for fixed y and that a minimum value of m must exist for some value of x. Since x 
appears on both sides of the above equation determining the horizon radius, a trivial relationship 
between m and x no longer occurs as it does for the ordinary D-dimensional Schwarzschild BH 
solution. In the 9,y ^ limit, corresponding to the usual commutative result, the upper limit of 
the integral defining the functions F^ becomes infinite and we arrive at the well-known standard 
result^ as then F„ ^ 1: 

m = c„x"+^ . (15) 

Note that if we had chosen a different form for the matter distribution representing the 
smeared point mass source only the set of functions Fn{z) would be changed but their general 
properties would be identical to those above. For example, if we had taken a smearing of the 
modified Lorentzian form, p ~ (r^ + 0'2^-{n+4)/2^ (with the parameter 6' not necessarily being the 
same as 9) then the corresponding functions, which we'll call G„, would also vanish as z — > in 
a power law manner and go — > 1 as z ^ 00 in a monotonic fashion. For example, one obtains 
Gq{z) = ;|(tan^-'^ z — Yqrp-) which we observe has these same limiting properties and is quite similar 
to the Fn qualitatively. More generally we find that the G„, are given by the integrals 
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These basic properties of the F„ (and Gn) capture the essential aspects of the NC physics. Other 

possible smearings, such as a straightforward exponential, would also lead to quite similar results.^ 
We will have more to say about this below but note that many of the specific expressions that we 
will obtain remain applicable if we take a different form for the smeared mass distribution. 

It is interesting to inquire how the NC value of x, for the moment explicitly denoted as 
xnc, compares with the usual D-dimensional result. Clearly for any fixed value of m the ratio 
xnc/x can be expressed solely via the functions Fn and is thus only dependent upon the ratio 
z = x/y and n; the result of this calculation is shown in Fig. ^ In examining these results, as well 
as those in the following figures, it is important to remember than when n = 0, = Mpi; for 
all larger values of n, ~ 1 TeV. In this figure we see that for large values of z <; 3 we recover 

the commutative result for all n since the NC scale is far smaller than the horizon size in this case. 
However when the two scales are comparable or when the horizon shrinks inside the NC scale the 
value of X is greatly reduced for fixed m. Of course, we really should not trust the details of our 
modeling of the NC effects when z is extremely small, i.e., when x is very much less than y. It is 
also in this very region where most of the differences between, e.g., the Gaussian and Lorentzian 
forms of the smeared mass distribution would be expected to begin to appear. 

More generally, we can now calculate m as a function of x for fixed y. Instead of a mono- 
tonically increasing function ~ we find that m[x) is now a function with a single minimum 

and which grows large as either x — > or ^ oo. For small x we find the scaling m ~ 
while the usual commutative behavior is obtained at large x, ~ x""^^. Fig. |21 gives a first rough 
indication of this behavior. The existence of a minima has several implications: (i) a minimum 
value of m implies that there is a physical mass threshold below which BH will not form, (ii) The 
inverse function, x{m), is double valued indicating the existence of two possible horizons for a fixed 

value of the BH mass; this is a potentiality first pointed out by NSS and something we will return 

^For the case of an exponential, the corresponding functions, Qn[z), are found to be related to the Fn above as 
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Figure 1: Horizon size in the NC scenario compared to the commutative result as a function of 
z = x/y for fixed values of m. On the left-hand side of the figure, from bottom to top, the curves 
correspond to n = 1 to 7. 

to below. In Fig. |21we see several additional features: first, as y increases for fixed n so does m 
except where x is large and we are thus residing in the commutative limit. Also for large y we see 
that m increases with n as it usually does in the commuting case. This is not too surprising given 
the small x scaling behavior of m above. Secondly, we see that for large y the values m are always 
large and the position of the mass minimum moves out to ever larger values of x as y increases. 
For example, if y = n = 1 then m ^ 400 and lighter BH do not form. Such large mass values are 
far beyond the range accessible to the LHC if we assume ~ 1 TeV, i.e., the interesting range 
roughly being 1 ^ m ^ 10 or so. ^ 

Since we are interested here in BH that can be created at colliders we will restrict our 
attention to smaller y values. In fact a short calculation shows that we need 0.05 ^ y ^ 0.2 in 

order to get into the LHC accessible mass region 1 ^ m ^ 10. To demonstrate this we must find 

the minimum value of the BH function of y for various values of n. This can be 

^Note that we do not expect BH to form for masses much less than Af, so that it is reasonable to believe that 
m > 1. 
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done in a two-step procedure: first we find where in x the mass minimum occurs for fixed values 
of y, i.e., where dm/dx = 0. Calculating this derivative we see that the minimum can be obtained 
by solving the equation (using for convenience the variable q = x^/4y^ introduced above): 



2g(n+3)/2g-<7 

(^ + i)r(^) 



^n(9)-7-^^^TT;7TIT3V=0' (17) 



which has a single, non-zero root qo{n). For any y this tells us the value of the horizon radius where 
the minimum mass occurs, Xmin = '^y\/ Qo{n), which we can now use to obtain rrimin employing 
the equations above. The result of this calculation is shown in Fig. |21 Here we see that for n in 
the range 1 to 7 the relevant values of y are rather narrow, not differing from y ~0.1 by more than 
about a factor of 2. 

It is important to recall that for the ordinary commutative L>-dimensional BH solution 
both rrimin and Xmin are algebraically zero. However, since we believe that rrimin ^ 1 is required 
to produce a BH, stronger conditions are usually imposed. In our case the results shown in Fig. |31 
represent the algebraic lower bound on m which certainly — > as — > 0. Physically, we might 
crudely imagine that rrimin ~ ^^^.^i^^^min) where rn'^-^ is the result shown in Fig. |21 We note 
that having a finite rrimin implies a BH production threshold while a finite Xmin implies a minimum 
BH production cross section, asH — '^^'^min/^* ^ colliders such as the LHC as is shown in Fig. ^ 
Here we see that far above threshold the BH production cross section scales like ~ m2/("+i) as 
would be expected in the commutative theory. However, for lighter BH this cross section falls 
significantly below this simple scaling rule and becomes quite small in the neighborhood of rrimim 
almost but not quite vanishing. 

We can now focus on this region of small y; for simplicity in what follows we will generally 
concentrate our results on the case y = 0.1. Fig. Elshows m{x) for y = 0.1 — 0.2 for n in the range 
to 7. Note that generally larger n leads to smaller horizon radii for fixed y and already at y = 0.2 
we see explicitly that BH are too massive to be produced at the LHC as expected from the above 
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Figure 2: BH mass as a function of the horizon size for n = O(soUd), n = 1 (dashed) and n = 
2(dotted). The upper (middle, lower) set of curves correspond to y = 10(1,0.1). 
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Figure 3: Minimum BH mass as a function of y showing the (very) approximate mass range 
accessible to the LHC between the dashed lines. On the left-hand side from top to bottom the 
curves correspond to n = 1,3,5 and 7, respectively. The allowed parameter range is above and to 
left of each curve. 
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Figure 4: NC BH production subprocess cross section as a function of m for y = 0.1. From right 
to left at the bottom of the figure the curves correspond to n=l to 7. 



analysis if ~ 1 TeV. One also sees that at x ~ 1 the asymptotic behavior, m 
already begun to set in for all n. 



,.n+l 



has 



As an example of the insensitivity of these results to our Gaussian parameterization of the 
smearing due to NC effects, let us briefly considered the modified Lorentzian form mentioned above. 

Setting y' = M^^/W = 0.1 (which is not necessarily the same as y = 0.1), we again evaluate m as 
a function of x. The result of this calculation is shown in Fig. El which we should compare with 
the top panel of Fig. El that yields essentially the same result in the asymptotic large x ^ 0.5(as 
it should since this is the commutative limit). In both cases a minimum mass occurs at relatively 
small X with somewhat similar values of m. rrimin decreases in both cases as n is increased and 
the corresponding value of Xmin also decreases as n is increased. The rrimin values are seen to be 
quite comparable in the two cases. The greatest difference in the two results is seen to occur in 
the region of x below Xmin where there is the most sensitivity to NC effects and the detailed shape 
of the BH mass distribution. This is just what we would have expected; in the region where NC 
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Figure 5: Same as in Fig. [^J but now for y = O.lO(top), 0.1 5 (middle) and 0.20(bottom) with the 
dotted curve being for n = 0. From top to bottom on the left-hand side of the figure the sohd 
curves are for n=l to 7. 



effects just begin to be felt the detailed nature of the peaked mass distribution is not actually being 
probed. What is really being probed in this parameter range is the fact that there is a peaked mass 
distribution instead of a 5-function source, i.e., the BH has a form factor due its finite size, and 
not the details of its shape. Only at smaller values of the radii (relative to the values of y or y') do 
the differences in the details of the mass distribution become important. In fact, we can if we wish 
tune our chosen value of y' to make these two sets of curves even more alike. Since the majority 
of the effects we will discuss are sensitive only to physics with x > Xmin, this short analysis shows 
that the general features of the results that we obtain below are not particularly sensitive to how 
the NC smearing of the mass distribution is performed. 




top to bottom in the middle of the figure the curves correspond to n = 0, 2, 4 and 6, respectively. 

In order to understand the possibility of the formation of two (or no) horizons, we follow NSS 
and consider the metric tensor component 5^00 as a function of the dimensionless radius co-ordinate 
M*r as shown in Fig. [3 Here we will assume that m = 5, a typical value which is kinematically 
accessible at LHC, for demonstration purposes. Recall that horizons occur when (700 = 0. With 
y = 0.1 all of the curves pass through 500 = twice corresponding to two horizons, one on either 
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side of Xmin- This explains why x{m) is double valued in Fig. [51 i.e., the two solutions correspond 
to the two radii where goo vanishes. For n = 0, the case studied by NSS, these two horizons are 
rather close in radius but this separation grows significantly as n increases. When y = 0.15, we 
see that for n = — 4 no horizon form as ^oo ^ 0.13. This corresponds to the result observed in 
Fig. El for y = 0.15 where we see that for this range of n the value m = 5 is not allowed. However, 
for n > 5, we again obtain two horizons; clearly a tuning of parameters will allow the two horizons 
to converge to the case of a single degenerate horizon at Xmin as found by NSS. For both values 
of y we note that as M*r — > the metric is no longer singular as in the pure Schwarzschild case 
as was noted by NSS in 4-d (as we are inside a well-behaved mass distribution), independently 
of the existence of any horizons. This is further confirmed by constructing the Ricci curvature 
invariant, R, as can be easily done from the Einstein equations above; in fact, we find that as 
X ^ 0, ~ -y-("+3) Furthermore, and more explicitly, apart from an overall numerical factor, 
R ~ my^("+^) [n + 4 — iMi^]e-(^^*^)^/%^ go that R is seen to vanish as M*r ^ oo as expected and 

undergoes a change of sign in the region, e.g., M^,r ~ 0.2 — 0.5 for y = 0.1 independently of m and 
only weakly dependent on the value of n. As we will now see only the 'outer' horizon, i.e., the one 
with X > Xmin is actually relevant to us. 

Our next step is to determine the thermodynamic behavior of these NC BH; to do this we 
first must calculate the Hawking temperature of the BH. This can be done in the usual manner by 
remembering that 



i.e., the temperature is essentially the r (radial) co-ordinate derivative of the metric evaluated at 
the horizon radius. Defining for convenience the dimensionless temperature, T = Th/M^, we obtain 
from the above form of the metric 



1 dd 



r=x/M:t 1 



(18) 



47r dr 



T = 



n + l 



1 - 



2g(n+3)/2g-g 



(19) 



47rx 



Fn{q){n + l)T{^) 
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Figure 7: goo as a function of M^r assuming m = 5 for either y = O.lO(top) or 0.15(bottom). The 
dotted curve corresponds to n = while from top to bottom on the left-hand side of the figure the 
solid curves are for n=l to 7. The dashed line corresponds to goo = 0. 
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Note that as expected T returns to the usual commutative result in the q ^ oo limit, i.e., the 
quantity in the large square bracket above ^ 1 in this limit. It is instructive to compare the 
temperature we obtain in both the NC and commutative cases for various values of the parameters; 
the ratio of these quantities is shown in Fig. |HJ This ratio is seen to be near unity for large z as one 
would expect but it decreases rapidly as z approaches the ~2-3 range from above. The temperature 
is also seen to vanish at the same z value where the BH mass is minimized. For smaller z, T becomes 
negative (which is where the second horizon occurs) and thus we enter a region that might usually 
be considered unphysical. If we had instead chosen the Lorentzian smearing these results would be 
quite similar qualitatively. 




1 2 3 4 5 6 7 

Z 



Figure 8: The ratio of the BH NC temperature to that obtained in the commutative limit where 
there is only a point matter source. The dotted curve corresponds to n = whereas, from bottom 
to top on the left-hand side of the figure, the solid lines correspond to n=l to 7. 

Fig. ini shows the actual NC temperature as a function of x for y = 0.1 and different n 
values. Here we see that T is quite close to its commutative value for x near unity, goes through 
a maximum as x decreases and then falls to zero at the BH mass minimum point. From these 
results and Fig. [S] above we can now trace the history of the entire semiclassical BH evaporation 
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process. " Consider a BH formed in a suitable parameter space region with moderately large values 
of m ~ 8 — 10 (and, hence, with large x). For such BH their Schwarzschild radii are too large to 
feel the effects of the NC scale in the formation process since x » y. As in the usual commuting 
picture, when the BH emits Hawking radiation it loses mass and gets hotter and thus radiates even 
more quickly. As the BH shrinks it begins to feel the NC effects and the temperature reaches a 
maximum in the mass region m ~ 1 — 7, the specific number depending on the value of n. As 
the BH continues to lose mass its temperature now decreases so that it radiates ever more slowly. 
Finally, as m approaches rrimin the semiclassical radiation emission processes ceases since T ^ 
has been reached leaving a classically stable remnant. Though sounding somewhat unusual such a 
possibility has been discussed in the literature for a number of alternative BH scenarios which go 
beyond the basic picture presented by General Relativity based on just the EH action pTl I13j . ** 
Whether quantum effects destabilize such a relic is not known. 

It is easy to convince oneself that a classically stable remnant is the natural outcome of this 
scenario. In the usual treatment of BH decays, the mass loss rate (assuming a perfect radiator) is 
given by 



with Hrf being a positive numerical constant. For decays dominantly to bulk(brane) fields we have 
d = n(0). In either case, clearly the lifetime of the BH is then given by 



with rriinitiai being the original BH mass. We recall from above that while x{m) is double valued 
it is well-behaved and never vanishes. However, on the otherhand, we also know that T(m) — > 

as 777. > Ttlmin- Thus for all n (and any d), Ibh will be driven to infinity due to the presence of a 

"Before beginning this discussion we must recall the usual argument discussed above that if rrimin lies below ~ 1 
then no BH will form; we will ignore this prejudice in the following discussion treating rrimin as the true minimum 
BH formation mass as derived from the theory itself. 

**These include models with higher curvature invariants in the action as well as those where a minimum length 
scale exists or where the Newton constant is taken to be a running parameter. 



dm 
~dt 



(20) 




(21) 
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singular denominator in the integrand implying a stable relic. (Of course, an initial very massive 
BH will radiate down to a mass very close to rrimin quite quickly.) This must happen in any model 
that predicts T ^ at finite x due to a corresponding singularity. For the commuting case, since 
T ~ 1/x, the integrand is never singular so that the BH lifetime remains finite. 

The lower panel in Fig. [^summarizes this discussion where we see that as m decreases from 
a large value the temperature increases, reaches a maximum value and then falls to zero at Tn^iyi 
leaving a classically stable remnant. 

This unusual temperature behavior in the NC case can also be studied more fully by ex- 
amining the BH heat capacity/specific heat, C. In the commutative case, C is always finite (away 
from X = 0) and negative since the BH gets hotter as it loses mass. Let we define 



C 



dm dm I dT 
'df ~~dx\~dx 



-1 



(22) 



which we can explicitly write as 



C = -47rc„x"+2F„((7)-i 



2H„iq) 



2Hn(q) , 4qH„{q) 
(n+1) (n+1) 



(n+3) 
2q 



Hn(q) 
1 



(23) 



where for convenience we have defined the set of auxiliary functions 



Hn{q) 



^(n+3)/2g-g 



(24) 



Using this expression, Fig. shows the results for the NC BH heat capacity as a function of 
both X and m for our standard choice of y = 0.1. We first see that C remains negative at large 
X and asymptotes to its commutative value, C — — Aircn x""*"^, as it should. Further we note that 
at X = Xmin (or at m = mmin) we find that C ^ as was expected. BH with mass mmin are 
no longer capable of mass loss since they have zero temperature. Between these two regimes the 
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Figure 9: The NC BH temperature as a function of x(top panel) and m (lower panel) for fixed 
y = 0.1. The dotted curve corresponds to case of n = whereas, from bottom to top on the 
right-hand side of the figure, the solid lines correspond to n=l to 7. 
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behavior of C is quite interesting. Consider C as a function of large and decreasing x. C at first 
decreases in magnitude as it does in the commutative case. However as we know from above, 
for some n-dependent x value, T reaches a maximum and then decreases. This implies that the 
magnitude of C then increases and becomes singular for this particular x value. For lower x the 
sign of C changes as now ^ > and then approaches zero at rrimin- This is even more obvious 
when we consider C as a function of m. For large m the BH radiates as in the commutative case as 
C decreases in magnitude as m is reduced. However, at some point the NC effects turn on and — C 
begins to increase, becoming singular where ^ = 0, the location of temperature the maximum. 
For smaller masses, as the BH radiates and gets lighter the temperature decreases so that we are 
now in a region of positive heat capacity. As m decreases further to rrimim C becomes zero for the 
remnant. 

Our next step in examining the thermodynamics of NC BH is to consider the value of the 
entropy which is defined via 

/dm yn+i 
^^^-'9^=4-c. / dv ^ , (25) 

where here we have made the natural choice that the entropy vanish at Xmin where the BH mass m 
is minimized for fixed values of y and n. In the commutative limit where 0,y ^ then Xmin — * 
and we recover the usual lower limit of integration. Fig. ^2 shows the values of the entropy for 
various n as a functions of x or m for our canonical choice of y = 0.1. We again observe that the 

commutative power law behavior, S = 47rc„ ^^q^gj > is recovered in the large x ^ 1 limit as expected. 
It is more interesting to consider S as a function of m; as we see from the Figure, while the entropy 
scales as ~ m^^^'^^^^'^^^^ for large m it rapidly falls below this scaling law to zero as m approaches 
nT-min from above. 

Finally, we also examine the free energy of the NC BH which is simply given by combining 
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Figure 10: Negative of the NC BH heat capacity/specific heat as a function of x(top panel) and 
m(bottom panel) for fixed y = 0.1. The dotted curve corresponds to case of n = whereas, from 
bottom to top in the middle of the figure, the solid lines correspond to n=l to 7. In the bottom 
panel, the parameter region with negative temperatures has been removed. 
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Figure 11: The NC BH entropy as a function of a;(top panel) and m (bottom panel) for fixed y = 0.1. 
From left to right (right to left) at the bottom of the top(bottom) panel the the curves are for n = 1, 
3, 5 and 7, respectively. 
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the above thermodynamic quantities: 

F = m-TS. (26) 

As seen in Fig. 1121 it too returns to its commutative value, F = m/{n + 2), in the hmit of large 
a; ^ 1 or large m. However, for x = Xmm it is clear that F = m for all values of y and n since both 
T and (by definition) S vanish at this point. Immediately above x = Xmin (or rrimin), F decreases 
slightly, until it matches onto the asymptotic ~ x""^^ behavior. For some parameter values F can 
even become negative in this intermediate mass regime. Generally we see the unusual behavior that 
m = rrimin is not the location of a minimum in F as either a function of x or m. This minimum 
lies near the values of x and m where T is maximized. 

3 Discussion and Conclusions 

Space-time noncommutativity and extra dimensions are both well-motivated ideas within the string 
theory context and it is natural for them to make their appearance felt as one approaches the 
fundamental scale. In more than four dimensions it is possible for this scale to be not far from ~ 
TeV and thereby address the gauge hierarchy problem. If this mass parameter is indeed this low it 
is likely that black holes will be produced at the LHC in sufficient quantities that their properties 
will be well measured. The occurrence of NC also at a similar scale could lead to significant 
modifications in the anticipated properties of these BH. Since a complete NC theory of gravity 
does not yet exist it becomes necessary to model the NC effects within the commutative General 
Relativity framework. 

Nicolini, Smailagic and Spallucci presented a physically motivated model of this kind in 4-d, 
where the essential aspects were the Gaussian smearing of matter distributions on the NC scale 
and the continued applicability of the EH action. They then went on to examine NC effects on 
BH physics. In this paper we extended this NC BH analysis in several ways: (i) we generalized 
the NSS study to the case of extra dimensions with a fundamental scale in the TeV range so that 
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the associated BH can be produced with large cross sections and studied in detail at the LHC. 
While much of the BH behavior observed in extra dimensions was similar to that obtained in 4-d, 
some significant modifications to the previously obtained 4-d results were observed. However, there 
appears to be an overall dominance of NC effects over those that arise due to the existence of the 
extra dimensions, {ii) We demonstrated that the essential physics induced by NC smearing is not 
particularly sensitive to the nature of the smearing function. In particular we explicitly showed 
that Gaussian and Lorentzian smearing lead to essentially the same behavior for the expected 
modifications of the BH mass-radius relationship due to NC effects. {Hi) We extended the NSS 
analysis to include several other thermodynamic quantities which are of interest in the study of NC 
BH such as their entropy, heat capacity and free energy. 

Perhaps the most important qualitative influence of NC on BH physics was already ob- 
served in 4-d by NSS, i.e., the existence of a classically stable remnant whose mass and radius are 
completely fixed by the NC scale and the number of dimensions. Within the framework of extra 
dimensions, if the fundamental scale is not too large then BH and their remnants will be copiously 
produced at the LHC and studied in detail. The observation of NC effects in the properties of these 
BH can open a new window on the fundamental theory of gravity and space-time. 
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